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ABSTRACT 

High-resolution X-ray spectrometers onboard suborbital sounding rockets can search for dark matter 
candidates that produce X-ray lines, such as decaying keV-scale sterile neutrinos. Even with exposure 
times and effective areas far smaller than XMM-Newton and Chandra observations, high-resolution, 
wide held-of-view observations with sounding rockets have competitive sensitivity to decaying sterile 
neutrinos. We analyze a subset of the 2011 observation by the X-ray Quantum Calorimeter instrument 
centered on Galactic coordinates l = 165°, b = —5° with an effective exposure of 106 seconds, obtaining 
a limit on the sterile neutrino mixing angle of sin 2 26 < 7.2 x 10 -10 at 95% CL for a 7 keV neutrino. 

Better sensitivity at the level of sin 2 26 ~ 2.1 x 10 -11 at 95% CL for a 7 keV neutrino is achievable with 
future 300-second observations of the galactic center by the Micro-X instrument, providing a definitive 
test of the sterile neutrino interpretation of the reported 3.56 keV excess from galaxy clusters. 

Subject headings: dark matter — Galaxy: halo - line: identification — neutrinos — techniques: 
spectroscopic — X-rays: diffuse background 


1. INTRODUCTION 

A variety of dark matter models predict photon pro¬ 
duction via dark matter decay, annihilation, or de¬ 
excitation. Some of these models predict mono-energetic 
photons with energies in the 1-100 keV range, prompting 
recent searches for lines in existing X-ray data from the 
XMM-Newton and Chandra observatories. Due to the 
well-understood atomic physics in this energy range and 
the ability to check the morphology of a potential signal 
against expectations from galactic dark matter halos, X- 
ray lines could provide unambiguous evidence for some 
models of astrophysical dark matter. 

There has been heightened interest in dark matter 
searches in the X-ray band following claims of an uniden- 


vations o f the Milk y W ay Galactic Center (MW GC) by 
Boyarsky et al. (2014a) finds a formal 5.7fj excess at the 
expected energy! However, the complexity of the GC 
makes modeling the background difficult, and because of 
the instrumental resolution of the observation they can¬ 
not rule out the possibility of the excess coming from K 
XVIII emission. 

A vigorous se arch has ensued, with v arious repor ts of 


non-detections: R iemer- Sorensen (2014|) in th e GC, Jel- 
tema & Profumo ( |2014| in the GC and M31, [Anderson 


et al. ~POTf) m g alaxies and galaxy groups, and Maly¬ 


shev et al. (2014) in dwarf spheroidal galaxies. There 


has been some debate as to how to best fit the contin¬ 
uum and of what the allowed flux of astrophysical lines 


seen in both galaxy and galaxy clus- 
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line is present at the > 3er level in three separate subsam¬ 
ples of data from both the MOS and PN instruments, and 
they als o detect it in Chandra ob servations of the Perseus 
cluster. |Boyarsky et al.| (|2014b|) reported a > 3cr excess 
around 3.53 keV in their spectral fits of XMM-Newton 
observations of the Perseus cluster and the Andromeda 
galaxy. In both cases, the width of the measured ex¬ 
cess is determined by the XMM-Newton and Chandra 
instrument response. Analysis of XMM-Newton obser- 
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ing Suzaku but find its spatial distribution in tension 
with expectations for decaying dark matter, and further 
they do not detect the expected scaled emission in ei¬ 
ther the Com a, Virgo or Ophiuchus clusters. |Carlson 


et al. (2015) performed a morphological study of the 
continuum-subtracted excess emission at 3.5 keV in the 
GC and Perseus. They find the GC excess spatial distri¬ 
bution incompatible with the expected DM distribution, 
and strongly correlated with the morphology of atomic 
lines from Ar and Ca with energies between 3-4 keV. 
The Perseus emission is correlated most strongly with 
the cool core emission, confirming the tension presented 
in the Suzaku Perseus observation. 

This set of observations demonstrates the challenges in 
searching for X-ray lines with the current observatories. 
New instruments are needed to improve the sensitivit y of 
searches for X-ray line emission from dark matter. |Bo-| 
yarsky et al. (2007) studied the optimal characteristics of 
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a mission dedicated to searches of diffuse line emission in 
the X-ray regime. They pointed out that the main deter¬ 
minants of instrument sensitivity are grasp = Aeftflpov 
(effective area x field of view, also referred to as etendue) 
and ene rgy resolution A E/E. As a “prototype” demon¬ 


stration, 


Boyarsky et al. (20071 calculated the limits on 
the sterile neutrino mixing angle sin 2 29 for data from the 
third flight of the X-ray Quantum Calorimeter ( XQC) 
sounding rocket payload (McCammon et al. 20021 for 
sterile neutrino masses between 0.4-2 keV. 

Existing X-ray telescopes tend to have comparatively 
small fields of view (e.g. a few tens of arcminutes for 
instruments on XMM-Newton and Chandra) with insuf¬ 
ficient energy resolution to resolve closely spaced weak 
spectral lines (e.g. ~ 100 eV FWHM at 2 keV for the 
EPIC camera on XMM-Newton) . Because of our location 
within the dark matter halo of the MW, the sterile neu¬ 
trino decay is an all-sky signal , so sensitivity can be im¬ 
proved by increasing the FOV. Discrimination of a signal 
against atomic lines can also be significantly improved 
using the superior energy resolution available with X- 
ray microcalorimeters. This combination of large-FOV 
with high spectral resolution is achieved in existing mi- 
crocalorinreter payloads on sounding rockets. Although 
the exposure from a typical sounding rocket flight is less 
than 300 s, the sensitivity of these short observations can 
be competitive with deep XMM-Newton observations of 
the GC. The upcomi ng SXS microcalorimet er instrument 
onboard ASTRO-H (|Takahashi et al.|2014 ) will have ex¬ 
cellent < 7 eV resolution, but its narrow 3’x3’ FOV 
limits its sensitivity to the all-sky signal expected from 
sterile neutrino decay in the MW. Wide-FOV sounding 
rocket observations are therefore complementary to the 
deep (~1 Msec) observations of the cores of galaxy clus¬ 
ters, gala xies, and dwarf s p hero idals that ASTRO-H will 
perform ( [Kitayama et al. 2014). 

In this paper we set limits on decaying sterile neutrino 
dark matter using a new dataset from the XQC sound¬ 
ing rocket in order to demonstrate the reach and analysis 
of large-FOV observations, and we discuss the sensitiv¬ 
ity and optimization of future observations using new 
instruments, such as the Micro-X detector. Section [2] 
discusses the sterile neutrino signal and estimates the 
signal and background of large-FOV observations for a 
putative signal. In Section [3| we describe the XQC in¬ 
strument, present an analysis of data taken during the 
5th flight of XQC, and place limits on the sterile neu¬ 
trino mixing angle sin 2 29 for sterile neutrino masses of 
4-10 keV. Section [4] estimates the sensitivity of observa¬ 
tions with the future Micro-X payload by constructing 
a detailed background model based on ASCA, ROSAT, 
and Suzaku observations and analyzing mock data sets. 
Section [5] discusses the implications of our flux limits on 
sterile neutrino dark matter. 


2. KEV DARK MATTER WITH ROCKETS 
2.1. Dark Matter Interpretations of X-ray Lines 


Well-motivated dark matter models can produce X-ray 
lines through decay, de-excitation, or annihilation. Per¬ 
haps the best -known scenario is tha t of keV-mass ster- 


ile neutrinos I 
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LJJ2006), although a large number ot mod- 


els have been proposed following the observations of the 


3.56 keV line. Such models include axions, axinos, excit¬ 
ing dark matter, gravitino^ mo duli, and WI MPs, a mong 
others (see discussion in Jeltema & Profumo [2014). 


Sterile neutrinos and other models that produce pho¬ 
tons by particle decay predict a flux per volume element 
that scales as the dark matter density ( p ). Oth er mo d- 


2007 

2014 

Berlin et al. 

2015 

), ti 


element that scales as the dark matter density squared 
(p 2 ). More complex scenarios, involving eXciting dark 
matter with a p rimordial population in an excited state 


(Finkbeiner & Weiner 2014), can furthermore produce 
fluxes that scale as p“ with 1 < a < 2. A non-linear 
scaling implies a much smaller flux from lower density 
systems like dwarf spheroidal galaxies and could ease the 
tension with non-observations of the 3.5 keV line in these 


systems (Malyshev et al.| 2014). In this paper we will use 

~D 


sterile neutrinos as our benchmark model, although we 
present our results as a line flux limit from a particular 
target, which can be translated into a constraint or signal 
in any of these models. 

Sterile neutrinos with masses in the ^ 1-100 keV range 
may contribute to the dark matter relic density if they 
are produced in the early universe. Two well-studied 
production mech anisms a re non-resonant os cillation of 
active neutrinos (Dodelson & Widrow 1 994]) an d reso - 
nant oscillation via the MS W effect ( Shi fc Fuller||1999 1. 
The existence of sterile neutrinos is additionally moti¬ 
vated by neutrino oscillation data, which could be ex¬ 
plained by adding sterile right-handed neutrinos to the 
standard model, such as in the jzMSM scenario (Asaka & 


Shaposhnikov 20051. Although they must possess cosmo- 
logical lifetimes in order to contribute to the dark matter 
relic density, sterile neutrinos may decay to a photon and 
active neutrino via a loop-suppressed process mediated 
by oscillation between the acti ve and sterile states. Th e 
rate for this process is given by (Pal fc Wolfenstein|1982 ) 
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where m s is the sterile neutrino mass and 9 is the mixing 
angle between the active and sterile states. 

Limits on sin 2 29 depend on both the observed flux and 
the fraction of dark matter comprised by sterile neutri¬ 
nos. For simplicity, limits are typically quoted under the 
assumption that sterile neutrinos comprise all of the dark 
matter. The X-ray fluxes reported by the claimed detec¬ 
tions discussed above in stacked galaxy clusters, M31, 
Perseus, and the MW GC correspond roughly to sin 2 29 
of 10" 11 to 10~ 10 . 


2.2. Dark Matter Signal for Large FOV Observations 

The flux expected from decay of sterile neutrinos in 
the MW halo is proportional to the integral of the DM 
density along the line-of-sight and over the field of view 

F=-^-l [ p(r(£, ip)) dl dfl, (3) 

m s 4t r J FOV J 0 
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distance from galactic center [kpc] 


Figure 1. Example of common DM halo profiles using NFW 
dNest i Sz Saluccil20131) (dotted), Ein asto ( Bernal & Palomares-Ruiz| 
[Mf| ~ (dashed), and cored Burkart (|Nesti Sc Salucc‘i| 2 U 1 2 >[: (sotid) 
parameterizations. 


where in the integral of the dark matter profile density 
p{r), the parameter l is the distance along the line of 
sight, r is the distance from the GC, and the angular 
integral is taken over the field of view of the instrument. 
The distance from the GC is related to the line-of-sight 
distance by 

r(£,ip) = \/£ 2 + d 2 — 2 Id cos ip, (4) 

where d is the distance of the earth from the GC and ip 
is the opening angle from the GC. 

In order to obtain a large FOV, the sounding rocket 
observations considered in this paper do not use an X- 
ray optic. The detector observes a field determined by 
an optical stop and has no imaging capability. The ef¬ 
fective area is that of the detector itself, on the order 
of 1 cm 2 . Furthermore, sounding rocket flights observe 
for a few hundred seconds per flight. In comparison, 
XMM-Newton has made observations on the order of a 
megasecond, with an effective area at 3.5 keV of around 
200 cm 2 for each MOS detector. 

In order to compare between different FOV observa¬ 
tions, a DM halo must be assumed, and we show several 
representative profiles in Figure [l] In Figure [2] we show 
in black the ratio of the expected rate from sterile neu¬ 
trino decay in the central 14’ radius of the GC (XMM- 
Newton’s FOV) to a larger FOV also centered on the GC. 
In gray we show the ratio between the XMM-Newton GC 
observation and a different field near the MW anti-center 
at Galactic coordinates l = 165°, b = —5°, observed by 
the 5th flight of the XQC (discussed in the next sec¬ 
tion). With a sufficiently large FOV, signal rates (in 
cts cm -2 s _1 ) 3 to 4 orders of magnitude larger than the 
XMM-Newton GC observation are attainable. The back¬ 
ground (in this case meaning all X-ray flux of non-DM 
origin) of observations that include the GC do not in¬ 
crease as quickly with FOV since the X-ray flux from the 
GC and the Galactic ridge are much stronger than the 
cosmic X-ray background (CXB) but extend to roughly 
±5° from the plane. This gives large FOV observations of 
the GC a better signal-to-noise ratio. Finally, the higher 
energy resolution of the microcalorimeter instruments in 
sounding rockets can cut the continuum background per 
energy bin by over an order of magnitude when compared 


Figure 2. Ratio of sterile neutrino decay signal in a large FOV 
to a 14’-radius FOV ( XMM-Newton ) around the galactic center as 
a function of the FOV half-opening angle. Black curves are for an 
observation centered at the GC, while gray curves are for an obser¬ 
vation centered on the XQC field of (h, l) = ( 165°.— 5°). Line style 
corresponds to diff eren t DM profiles: N FW d Nesti & : Salucci 2013| 
(dotted), Einasto (Bernal & Palomarcs-Ruiz 20 l^jl ( dashed), and 
cored Burkart (jNcsti & Salucci 2013) (solid) parameterizations. 
The FOVs for Micro-X and AQC are indicated by arrows on the 
horizontal axis. 


reference flux (in 14’ of GC) 

2.9 X 10 -5 cm -2 s -1 

scaled flux (in 20° of GC) 

6.1 x 10 _2 cm _2 s _1 

effective area at 3.55 keV 

1 cm 2 

exposure time 

300 s 

resolution (FWHM) 

3 eV 

signal events (in 20° of GC) 

18.2 

bg. rate at 3.55 keV (see (4.1 

4.5 cm -2 s _1 keV -1 

bg. events in signal window 

6.7 

median signal significance 

5.6 a 


Table 1 

Basic signal and rates expected for a hypothetical observation of 
the GC using a microcalorimeter on a sounding rocke t, assuming 
a fiducial signal flux from|Boyarsky et al.| |2014a^. 


to the CCD energy resolution of XMM-Newton. 

To get a feel for the potential sensitivity of these obser¬ 
vations, consider a hypothetical instrument with a 1 cm 2 
effective area, 3 eV FWHM resolution, and 20° radius 
FOV at 3.55 keV, which observ es the GC for 300 s (see 
Table fTj) . If the flux reported by Boyarsky et al.| (2014a I 
over a 14’ radius FOV is due to decaying sterile neutrino 
dark matter with the NFW profile of Figure [l] then our 
hypothetical 20° field would expect a scaled signal flux of 
6.1 x 10 _2 cm _2 s _1 , 2000 times higher than the XMM- 
Newton observation. A 300 s observation would measure 
18.2 total events in the X-ray line. A t 3.5 keV, the back¬ 
ground model described in Section |4.1| predicts a flux 
of 4.5 cts cm -2 s _1 keV -1 . The background in a 5.1 eV 
(±2 <je) window would be 6.7 events. In spite of the small 
statistics, the median significance of the putative signal 
above the continuum background would be 5.6cr in this 
short observation. 

For comparison, in the (Boyarsky et al. (2014a) analy¬ 
sis of 1.4 Ms of XMM-Newton data we estimate around 
7,500 signal counts in the claimed 3.54 keV line in each 
MOS detector. In that same resolution element, there are 
upwards of 500,000 background counts. With a signal- 
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Figure 3. All-sky X-ray map from the MAXI/GSC instru ment 
onboard the International Space Station (|Mihara et al.|2014|). The 
image is a negative made from a color rendition with the loliowing 
energy band definitions: red, 2-4 keV, green, 4-10 keV, and blue, 
10-20 keV. The XQC field analyzed in this section centered on 
l = 165°, b = —5° is delineated by the solid line. The dotted line 
is a Micro-X field centered on l = 162°, 6 = 7°, chosen to lie inside 
the XQC field and evade the Crab pulsar. The dashed line is the 
Micro-X GC field and the long-dashed line is the Micro-X off-plane 
field centered on l = 0°,6 = —32°, both discussed in Sectional 

to-noise ratio of 0.015, the authors use the XMM-Newton 
high statistics measurement to detect such a small signal 
at high formal significance (5.7a), but doing so depends 
on an accurate model of their background and minimal 
systematic errors. 


3. ANALYSIS OF XQC DATA 

Having laid out the basic strategy, we now focus on 
existing data from X-Ray Quantum Calorimeter. The 
XQC payload is a mat ure f light system wit h 6 flights be¬ 
tween 1995 and 2014 (Crowder et al. 2012). The XQC 
spectrometer is an array of 36 microcaiorimeters using 
ion-implanted semiconductor thermistors each coupled 
to a 2 mm x 2 mm x 0.96 /jin HgTe absorber on a 
14 /un-thick Si substrate, with total area of 1.44 cm 2 . 
The energy resolution below 1 keV is 11 eV FWHM, 
although due to position dependence the resolution de¬ 
grades to 23 eV FWHM at 3.3 keV. The microcalorime¬ 
ter array is mounted inside a cryogenic system which 
uses pumped He as a 1.5 K bath for an Adiabatic De¬ 
magnetization Refrigerator (ADR), which is coupled to 
the detector assembly and cools it to ~ 50 mK temper¬ 
atures. To survive the launch vibrations while cold, the 
cryogenic system is suspended with vibration insulators 
from the skin of the rocket, and the resonant frequencies 
of the system are designed to minimize coupling of skin 
vibrations to the detectors during launch. The FOV of 
XQC is 1 sr, subtending a 32.3° radius in the sky. 

Data from an observation centered at the Galactic co¬ 
ordinates of l — 90°, b — 60 ° during the 3rd f l ight o f 
XQC were first presen ted by McCammon et al. (2002). 
Boyarsky et al.| (|2007|) used this data to constrain the 
decay of sterile neutrino dark matter, and their results 
are shown in Figure |12| Their analysis did not perform 
background subtraction and was limited to data below 
~ 1 keV. 

We perform a new analysis that develops a background 
model for the data between 2.0 keV and 5.0 keV, and 
then uses the data and background model to constrain 
the flux of an unidentified line in this interval. The use 
of background subtraction and higher-energy data from 
a more recent flight of XQC are the main improvements 


over Boyarsky et al. (2007). We analyze a partial data 
set from the fifth flight of the XQC rocket, which flew 
2011 November 06 at 08:00 UT as flight 36.364UH from 
the White Sands Missile Range. It obtained about five 
minutes of on-target data at altitudes above 160 km. 
The field of view was centered at the Galactic coordi¬ 
nates of l = 165°, b = —5° (shown in Figure |3|, close to 
the galactic anti-center and including the Crab nebula. 
A total of 200 s of on-target data was analyzed on 29 
functional pixels. After a very conservative quality cut 
to remove pixels and time periods with unstable event 
rates, 2551 pixel • s remain on 24 pixels, for an effective 
exposure of 106 s per pixel. Data from other XQC flights 
are also being reprocessed, and the combination of these 
data sets will increase the total exposure by a factor of 
a few in a future analysis. 

Figure [4] shows the XQC data above 2.0 keV. The spec¬ 
trum contains a power law continuum with strong lines 
at 3.31 keV and 3.59 keV from Ka and K/3 transitions 
of potassium, respectively. These lines arise from a 41 Ca 
source which provides continuous calibration during the 
flight, and is used to correct gain fluctuations. X-rays 
incident on XQC may be absorbed in either the HgTe ab¬ 
sorber or its Si substrate. The photons absorbed in the 
HgTe are efficiently thermalized, while those absorbed in 
the substrate experience an energy loss of about 15-20%. 
Potassium Ka and K fj events in the absorber’s Si sub¬ 
strate form the two broad peaks centered at 2.8 keV and 
3.0 keV, below the corresponding lines due to absorption 
in the HgTe. The relative intensity of the full-energy 
peak in the absorber and the second peak from events in 
the substrate is determined by the relative absorption ef¬ 
ficiencies for X-rays in the two detector elements, shown 
in Figure [5] XQC was optimized for studying the soft X- 
ray background in the 0.1-1 keV energy range, where the 
HgTe absorber has higher efficiency and almost all X-rays 
are absorbed before reaching the Si substrate. A detec¬ 
tor with thicker HgTe absorbers would be more suitable 
for the 2-5 keV region that we study here. The efficiency 
of the Si rises rapidly above 1.0 keV, and becomes com¬ 
parable to the HgTe efficiency above 3.5 keV. 

The overall strategy of the analysis is to perform an 
exclusion on the rate above background of an uniden¬ 
tified line centered at each energy between 2.0 keV to 
5.0 keV. To do this, we first fit a background model to 
the data incorporating the important spectral features 
described above. We then use the model to estimate the 
background in a sliding energy window, allowing us to set 
an upper limit on the expected flux from an unidentified 
line, as a function of energy. Modeling of the energy spec¬ 
trum becomes more complex at energies below 2.0 keV 
due to weak atomic lines from thermal emission. At en¬ 
ergies significantly above 5 keV the energy scale could 
become nonlinear and the detection efficiency drops due 
to saturation. For simplicity, we restrict to a conserva¬ 
tive energy window of 2.0-5.0 keV, although this range 
could likely be expanded in future analyses. 

The background model and its components are shown 
in Figure [4j The model incorporates the lines from the 
41 Ca source, a continuum from events in which photo- 
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Figure 4. Spectrum of XQC data overlaid with fitted total background model (solid blue). Dashed lines show background model 
components, consisting of a power law continuum from the diffuse X-ray background (long dashed cyan), a power law continuum from the 
Crab (dot-dashed red), cosmic rays (dotted purple), and lines from the 41 Ca calibration source onboard the instrument (short dashed green). 
The calibration source produces lines at 3.31 keV and 3.59 keV from Ka and K/3 transitions of potassium, while the two broad peaks at 
lower energies are due to Ka and K/3 X-rays that interact in the Si substrate of the HgTe absorbers and experience energy losses due to 
charge trapping. The flat continuum visible below 2.5 keV in the calibration spectrum is due to source events in which the photoelectron 
escapes the absorber. The bottom panel shows residuals between the data and the total background model, normalized by the error for 
each bin. 



Figure 5. Efficiency of XQC for detecting X-rays as a function 
of energy. Curves show the efficiency of absorption in the HgTe 
absorber (dashed orange), and the efficiency of absorption in the Si 
substrate (solid blue). Events absorbed in the Si substrate lose 15- 
20% of their energy due to charge trapping, but otherwise appear 
as good pulses (see text for discussion). 


The power law models are properly corrected for the sup¬ 
pressed energy measurement when X-rays interact in the 
Si substrate, as well as for the efficiency of X-ray de¬ 
tection shown in Figure [5] Fits are performed using 
the RooFit softwar e pack a ge based on the Minuit nu¬ 
merical minimizer (Verkerke & Kirkby 2003). We use 
the met hod of extended unbinned maximum likelihood 
(Barlow 1990), which is more stable than binned fits in 
a low-statistics setting. More details on the model and 
statistical methodology are provided in the Appendix. 

For each energy Eq, we construct a window [Eq — 
2a, Eq + 2a] of four standard deviations in energy resolu¬ 
tion which we use to set upper limits on any flux above 


the modeled background. The background model is fit 
to all data outside this signal window, and then extrap¬ 
olated into the window. The background model in the 
window is then integrated to obtain the background rate 
b with uncertainty erf, propagated from the uncertainty on 
the fit parameters. An upper limit is then set on the rate 
of signal above background in the window, for a Poisson 
process. Limits are s et using the profile likelihood test 

(2011), which incor- 


statistic described in Cowan et al. 


porates the uncertainty on the background rate in the 
window. The critical value of the test statistic for the 
desired confidence level is calculated exactly from Monte 
Carlo simulations rather than using the asymptotic dis¬ 
tribution of the test statistic. This is important at higher 
energies, where the low statistics cause the test statistic 
to differ from its asymptotic distribution. The up per lim¬ 


its are set^ using the RooStats software package (Moneta 


et al. 2011). Although the “sliding window” approach 


does not exploit the signal and background shape within 
the energy window, there is little loss of information be¬ 
cause of the very low statistics of a potential signal rela¬ 
tive to the slowly-varying background within the narrow 
window. 

Figure [6] shows the limits on the flux of an unidentified 
line as a function of the line energy, calculated using the 
limit-setting procedure described above with the back¬ 
ground model. At 3.53 keV, we set an upper limit on the 
flux of an unidentified line of 0. 17 cm~ 2 s -1 at 95% CL . 
The flux reported in the GC by Boyarsky et al.|( 2014a), 
when referred to the XQC field using fiducial L)iVl pro- 
files is listed in Table [2] The inset of Figure [6] also shows 
a strong downward fluctuation of the limit at 3.55 keV, 
consistent with the lower edge of the ±2cr band of ex- 
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Energy [keV] 


Figure 6. Upper limit at 95% CL on the flux of an unidentified line in the XQC spectrum (black line), as a function of the line energy. 
Bands show the ± I rr (green band) and ±2 <j (yello w band ) range of the exp ected limit from the background-only hypothesis. Inset shows 
the energy range around 3.5 keV. The flux from the|Boyarsky et al. (2014a| claim referred to the XQC field of view using the NFW profile 
from Figure[2]is shown with a red dot. 


Data 

DM profile 

Flux [cm 2 s - 1 ] 

XQC (this work) 

N/A 

< 17 x 10~ 2 (95% CL) 

Expected line flux in 
XQC field scaled from 

NFW 

Einasto 

Bur kart 

(2.03 ±0.35) x 10" 2 
(1.95 ± 0.34) x 10" 2 
(11.7 ±0.2) x 10“ 2 

Boyarsky et al. ( 2014a) 


Table 2 

Flux limits on a line at 3.53 keV from XQC data in this work, 
compared with the expected flux in the X QC fie ld obta i ned by 
scaling the galactic center observation of| Boyarsky et al. ([2014a), 
using the DM profiles of Figure Jii] 


pected limits. While this fluctuation may be purely ran¬ 
dom, it could also be caused by a very slight difference 
in the lower tail of the K/3 line between the flight data 
and the calibration template constructed using ground 
calibration data. Future flights of XQC would clearly 
benefit by choosing a calibration source with lines fur¬ 
ther from the signal region near 3.5 keV. Although the 
XQC data do not exclude the Boyarsky best-fit flux, the 
upper limit is close in the case of a cored DM profile, 
such as Bur kart. It is furthermore important to empha¬ 
size that XQC achieves this result with merely ~ 106 s 
of data on only a subset of its pixels. This underscores 
the value of combining additional datasets from other 
flights: additional statistics will both improve the limit 
and enable more robust background modeling. 


4. ESTIMATES FOR FUTURE OBSERVATIONS 

The XQC limits shown in the previous section are 
photon-limited; more observation time will result in bet¬ 
ter sensitivity. We also plan future observations in sev¬ 
eral fields including the GC to both increase sensitivity 
to potential undiscovered sterile neutrino lines and pro- 
vide a definitive test of the nature of the line found by 


Boyarsky et al. (2014a). These future observations will 
benefit from higher resolution microcalorimeter arrays, 



Figure 7. Total efficiency of Micro-X for X-ray detection as a 
function of energy. The total area of the detector is 1 cm 2 . 


such as those developed for the Micro-X rocket payload. 

The Micro-X payload is a new system based on the 
XQC design, with new detectors and readout to allow 
f or a larger array of higher-resolution microcalorimeters 
(Heine et al. 2014). Although Micro-X was designed 
to be used with a 2.1 m X-ray optic, for dark matter 
searches it will be reconfigured to fly without it, using an 
optical stop like XQC to set its field of view. The Micro- 
X detector consists of an array of 128 microcalorime- 
ters using Transition-Edge Sensor (TES) thermometers 
each coupled to a 0.6 mm x 0.6 mm BiAu absorber. 
For the dark matter flight, new absorbers with 0.9 mm 
per side and thickness of 3 /im (Bi) + 0.7 /im (Au) will 
be used, with total area of 1 cm 2 . The absorption effi¬ 
ciency of X-rays in the the detector is shown in Figure [7] 
The BiAu absorbers allow for better thermalization and 
should minimize position dependence, allowing Micro-X 
to retain its design resolution of 3 eV FWHM at 3.5 keV. 

In this section we estimate the sensitivity of a poten¬ 
tial GC observation with the Micro-X payload. A GC 
observation is the most direct comparison to |Boyarsky] 
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Figure 8. Expected X-ray background from a 20° radius region 
centered on the Galactic Center. The total background spectrum is 
shown in black. The expected emission from the brightest low mass 
X-ray binaries (as shown in Table |3[> is shown in green and labeled 
(1), the emission from the cosmicA-ray background is shown in 
blue and labeled (2), the thermal components from the galactic 
diffuse background emission are shown in green and labeled (3), 
and finally the purple line labeled (4) shows the contribution from 
ionized cold ISM neutral Fe, which is a combination of a powerlaw 
continuum and gaussians representing the Ka and K (3 iron line 
emission. The energy bins are 3 eV wide. 


et al. (2014a), and the higher energy resolution of the 
Micro-X instrument results in a lower background from 
continuum emission and better discrimination between 
unexpected lines and those coming from atomic transi¬ 
tions in the observed plasma. Micro-X is designed to 
be coupled to a mirror, so for this calculation we use a 
0.38 sr, 20° radius FOV as an estimate of the achievable 
FOV without the Micro-X optics. The Micro-X GC field 
is shown in Figure [3] A future redesign of the optical 
aperture of the cryostat could increase the FOV to 1 sr. 

4.1. Galactic Center Backgrounds 

The first step in estimating the sensitivity of a poten¬ 
tial MW GC observation with Micro-X is constructing 
a background model of the complex emission from the 
large FOV. For this we have used results from Suzaku 
observations of the Galactic Ridge (GR) and GC to es¬ 
timate the contribution from thermal diffuse emission 


(Ucliiyama et al. 2013), the cosmic X-ray background 
(CXB) compon ent from unresolved e xtragalactic sources 
as modeled by (Kushino et ah] "2002) from ASCA obser¬ 
vations, and the RUSAT All Sky Surv ey - Bright Source 
Catalogue (RASS-BSC, revision 1RXS, Voges et al.|1999|) 
to estimate the contribution from point sources in the 
field. The background model shown in Figures [8] and [9] 
is discussed below. 


4.1.1. Diffuse Background 

The spectral model for the diffuse emission expected 
from the field is constructed using the thermal compo¬ 
nents from observations of the GC and GR taken with 


Suzaku (Uchiyam a et al. 201 3), a s well as the isotropic 
CXB model from Kushino et al. (2002). The GC com¬ 


Figure 9. Background from a 20° radius region centered on the 
GC (as in FigureEl, in the 3.3-3.7 keV range. The total background 
spectrum is shown in black, and the components are as in Figure [8] 
Emission lines from ions Ar XVIII, S XVI, Cl XVI, Cl XVII, K XVIII, 
and Ar XVII are also shown. The abundanc e of Cl and K have 
been s et to solar values, using the tables from|Anders & Grevesse| 
(1989jl. The energy bins are 3 eV wide. 


ponent is defined as emission from a region with radius 
0.6°, centered at l,b = 0°,0°. The GR region extends 
from —5° to +5° in Galactic latitude and spans the full 
longitude range of the FOV (from 340° to 20°). We have 
defined the spatial extent of these regions using two- 
exponential model fit to the 2.3-8 keV intensity profi le 
around the GC presented in Table 2 of [Uchiyam a et al.| 
(20131. The integrated background spectrum is the sum 
of the emission from the GC, the GR, and the emission 
from latitudes \b\ > 4.9° to the edge of our FOV. 

Each of the Galactic components is a combination 
of emission from two thermal plasmas with different 
temperatures and metallicities, all absorbed through the 
appropriate column densities. The thermal plasmas are 
simulated in XSPEC (version 12.8.2 1) using the APEC 


model from AtomD B version 3.0.1 (Smith et al. 2001, 
Foster et al. 2012). We have limited our analysis to 


energies above 2.3 keV to avoid the contribution from 
additional (|&| > 5°) Galactic thermal components from 
regions extended beyond the GR region, since detailed 
spatial and spectral models of such extended diffuse 
thermal emission are currently not available. Neutral 
atoms from the cold ISM are ionized by either low 
energy cosmic rays or X-ray emission from external 
sources. This results in an additional diffuse X-ray 
component from the GC and GR, with a continuum 
modeled as a powerlaw (F = 2.13), and fluorescent 
K-shell line emission from neutral Fe (shown in purple 
in Figure |8j and Figure [9|. We used the normalization 
of the continuum and equ ivalent widths of the Fe Ka 
and K/3 lines presented in |Nobukawa et al.| (|2010|) and 
Uchiyama et al. (2013). However, we split the Fe Ka 
contribution into Fe Kal (at 6.403 keV) and Fe Ka2 


(at 6.390 keV), with 2:1 re lative intensities (|Bearden 
1967 |Kaastra & Mewe||l993]). 

lectral 


The CXB contribution is 


added as a powerlaw spectral component, given by 8.2 x 
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1 2 s 1 arcmin 2 keV 


10 _ ‘{E/1 keV) -1 ’ 4 photons cm - 
absorbed throu gh the appropriate co lumn density in 
each direction (Kushino et al. 2002 ). All spectral 
param eters are taken trom Table 3 of [Uchiyama et al.| 
(2013|), and the normalizations of all components are 
adjusted to match the nominal intensities predicted in 
that work for the spatial extents described above. The 
total thermal emission is shown in red in Figure [Sj and 
the integrated CXB emission in presented in blue. 


4.1.2. Background from Bright Sources 

The RASS-BSC includes 558 sources within 20° of the 
GC, with a total count rate in the ROSAT broad band 
(0.1-2.4 keV) of approximately 760 counts/s. The 12 
brightest sources within the region of interest, all low 
mass X-ray binaries (LMXBs), account for 80% (~ 613 
counts/s) of the total count rate from the resolved bright 
X-ray source population in the field, and we have listed 
them in Table [3] Their spectra are modeled using the 
emission and absorption parameters obtained through 
X-ray observations, using the references also included in 
Table [3j Note that the observed flux and spectral shape 
of these sources will depend on their state at the time of 
observation, so there is some inherent uncertainty in esti¬ 
mating this background component. Since the remaining 
portion of the X-ray bright source population is also dom¬ 
inated by LMXBs, the integrated spectrum from these 12 
brightest sources has been normalized in order for the to¬ 
tal flux in the 0.1-2.4 keV band to match the combined 
ROSAT count rate from all resolved bright sources in 
the region of interest. The integrated BSC spectrum is 
shown in green in Figure [8j 


4.1.3. Combined Sky Background 

Figure [8] shows the total background model for a 20° 
radius fiem centered on the GC, which includes the dif¬ 
fuse and point source emission discussed above. We show 
the total flux in black, the expected emission from bright 
sources (as shown in Table [31 in green, the CXB power- 
law in blue and thermal components from the Galactic 
Ridge and Galactic Center in red. Figure [9] shows the 
these different components but zooms in on the 3.3 keV 
to 3.7 keV region (line colors are consistent with those of 
Figure [8]). The continuum in this energy band is domi¬ 
nated by the emission from the BSC LMXBs, yet some 
thermal line emission is also significant. The relevant 
ions in this energy range are Ar xvm, S XVI, Cl XVI, Cl 
XVII, K XVIII, and Ar XVii. We u sed the abundances ta¬ 
ble of [Anders & Grevesse (1989) and set the abundance 
of Cl and K to those values. 


4.1.4. Instrumental Backgrounds 

In a realistic flight, Micro-X requires an on-board cali¬ 
bration source similar to the one that produces the potas¬ 
sium lines in Figure EJ The 41 Ca source used by XQC 
is not ideal for searching for a line at 3.5 keV because 
the KB line is at an energy similar to the signal. We 
consider an alternate calibration source consisting of an 
55 Fe source that produces fluorescence X-rays by illu¬ 
minating an NaCl wafer. A kapton filter could block 
Auger electrons from the NaCl and the ~ 1 keV X-rays 
from Na fluorescence, leaving only the Ka (2.62 keV) 
and K/3 (2.82 keV) lines from Cl fluorescence, as well as 



Figure 10. Mock data in the energy range of interest for the 
3.5 keV line, with (red) and without (black) a signal, with back¬ 
ground model (blue line) and signal model (dashed green line) 
overlaid. Note that the excellent spectral resolution of Micro-X 
provides significant separation of a signal line from nearby atomic 
lines. 


a small number of back-scattered 55 Fe X-rays. We sim¬ 
ulate th e energy spectr um obs erved by Micro-X using 
Geant4 (Agostinelli et al. 2003), and add it to the astro- 
physical X-ray spectrum, assuming a calibration source 
rate of 1 Hz / pixel. 

Cosmic rays will also produce a background in Micro- 
X. Most cosmic ray primaries are protons with energies 
around a few GeV. These act as minimum-ionizing par¬ 
ticles producing a broad continuum of energies in the 
detector, but peaking in the signal region around 3- 
4 keV. We simulate the cosmic ray energy spectrum using 
Geant4, and add this as a background component in the 
background spectrum. The total rate of cosmic rays is 
about 1 Hz in Micro-X, which is significantly less than 
the total rate expected from astrophysical sources. 

4.2. Signal and Sensitivity Estimates 

A mock observation o f the GC with and without the 
Boyarsky et al. (2014a I line is shown in Figure 10 A 
candidate line of this strength in the Micro-X obser¬ 
vation would also be stronger relative to atomic lines 
than in XMM-Newton and Chandra observations. In 
fact, since the background model for the GC is domi¬ 
nated by a blackbody continuum from low-mass X-ray 
binaries, no significant atomic lines are expected to be 
visible above the continuum in the energy range of inter¬ 
est. The strongest lines between 3.5 and 3.6 keV come 
from K XVIII and Cl XVII with expected counts of less 
than 1 event in the Micro-X observation, so if a line of 
this strength was detected at this energy it would be 
in strong tension with a “standard” astrophysical origin. 
As can be seen in Fig [lOj the putative sterile neutrino 
line would be bigger than the Ar XVII line at 3.683- 
3.685 keV (which is the sum of two emission lines), and 
this Ar line, in turn, is expected to be a factor of 30 (5) 
larger than the brightest Cl XVII (K XVIII) line. 

Using the background model and the methodology of 
our analysis of XQC data, we can estimate the sen¬ 
sitivity to an unidentified line over the entire energy 
range available in observations by Micro-X. We consider 
a 300 s measurement centered on the GC. While the 
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Table 3 

LMXBs included in the background model for the Micro-X observation. 


ROSAT Name 

Associated Name 

ROSAT Count Rate 
(counts/s) 

N h 

(10 22 cm 2 ) 

Fq .5 — 2 keV F2-IO keV 

(10 -9 ergs/s/cm 2 ) 

1RXS J173143.6-165736 

GX 9+9 

144.5 

0.15 

1.3448 

4.5127 

1RXS J182340.5-302137 

4U 1820-30 

127 

0.078 

1.3219 

5.188 

1RXS J170544.6-362527 

GX 349+02 

91.82 

0.7 

1.8785 

12.075 

1RXS J173858.1-442659 

4U 1735-44 

61.88 

0.185 

1.1257 

4.5544 

1RXS J175840.1-334828 

4U 1755-33 

31.08 

- 

- 

- 

1RXS J180132.3-203132 

GX 9+1 

30.92 

1.447 

0.71644 

19.997 

1RXS J173602.0-272541 

GS 1732-273 

29.79 

0.67 

0.82566 

1.2021 

1RXS J181601.2-140213 

GX 17+2 

24.82 

3.18 

0.6595 

15.843 

1RXS J170855.6-440653 

4U 1705-44 

20.37 

1.5 

0.5785 

6.4055 

1RXS J174755.8-263352 

GX 3+1 

19.2 

1.7 

0.42601 

4.4799 

1RXS J180108.7-250444 

GX 5-1 

17.78 

2.8 

0.6671 

41.512 

1RXS J173413.0-260527 

KS 1731-260 

14.09 

1.08 

0.40815 

2.9953 


Reference 



|Ng et al.| 

2010 



Costantini et 

al. (2012) 



Ng et al.j 

2010 

) 




Ng et al. 

2010 

) 


Ana; 

elini & W 

hite ( 

2C 

03) 

Iaria et al. 

(2005 


Yam 

auchi & Nakamura 

2004 


Cackett et al. 

(200< 



||Ng et al. 

20101) 


Piraino et a 

. 

(2012) 


jUeda et al. 

2005 



Narita et al. 

(2001) 
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Figure 11. Expected limit from an observation of a 20° field 
around the GC by Micro-X. Black line shows the median expected 
95% CL upper limit, while the green and yellow bands are the 
±lcr and ±2tr ranges of the expected upper limits. The expected 
upper limit rises at high energies because of the falling effic ienc y to 
detect x-rays. The red point is the flux of Boyarsky et al. |2014a), 
extrapolated to the Micro-X field of view using an .\ l ! W 7 profile. 
Gray bands overlay strong calibration source lines. 


Micro-X constraints are qualitatively similar to those of 
XQC, they are more stringent because of the significantly 
higher spectral resolution and exposure of Micro-X. 

Using the same analysis approach that we applied to 
XQC, we compute the expected upper limit on the flux 
of an unidentified line as a function of energy, under 
the background-only hypothesis. The resulting limit is 
shown in Figure [TI] Fluctuations of the limit at low en¬ 
ergies are due to the presence of atomic lines. At higher 
energies, large jumps in the limit are caused by the small 
number of background events in the signal region, while 
small fluctuations are due to finite statistics in the Monte 
Carlo simulation used to set the upper limit. 

5. STERILE NEUTRINO INTERPRETATION 

The flux limits obtained from XQC data and projected 
for Micro-X can be translated into constraints on models 
of dark matter. Although the literature contains a range 
of models that could produce an X-ray line, we consider a 
decaying sterile neutrino as a benchmark model because 
it has been extensively discussed as the source of the 


3.5 keV excess. Using equations ([2| and (|3j for the flux 
of a decaying sterile neutrino, and assuming an NFW 
profile, we translate the XQC flux limits of Figure [6] into 
limits on the sterile neutrino mass m s a nd m ixing angle 
sin 2 20, shown as the black line is Figure |l2[ 

For the Micro-X payload, we show three sensitivity 
projections, corresponding to the limits obtained for sim¬ 
ulated background-onlv_300 s observations of the three 
fields shown in Figure [3] The dotted line labeled (1) in 
Figure [12| shows the limit expected from a Micro-X ob¬ 
servation inside of the XQC field at l = 162° ,6 = 7°, 
chosen to avoid the Crab pulsar. For this projection, we 
take the background to be the diffuse cosmic X-ray back¬ 
ground measured by the XQC observation in Figure 0J 
scaled to the smaller Micro-X field of view, and witm 
out the Crab component. The dashed line labeled (2) in 
Figure [12| shows the limit expected from an observation 
of t he galactic center using the flux limit shown in Fig- 
ure [TT| Finally, the long-dashed line labeled (3) in Fig- 
ure |12| shows the limit of an off-Galactic-plane pointing at 
l = 0°, b = —32°. For this observation, the background is 
assumed to be only the cosmic X-ray background shown 
in blue in Figs. [8] and [9l reducing the continuum by a 
factor of ~9 (at 3.5 keV). The putative sterile neutrino 
signal is reduced by a factor of ~2 relative to the GC 
assuming an NFW profile, giving this field an improve¬ 
ment in signal-to-noise of 4.5 over the GC pointing. The 
limit, however, does not improve appreciably over the 
GC limit. This is due to the fact that the short 300 s ob¬ 
servations are photon starved, and the limits are driven 
by low-number statistics. Thus integrating over multiple 
flights will yield sensitivity improvements that will in¬ 
crease faster than sqrt(exposure). The off-Galactic-plane 
fields are particularly appealing since they have a much 
simpler background with fewer systematics than the GC. 
In the case of a positive signal, multiple pointings could 
be used to map out the expected DM profile. 

Although the XQC limit is not strong enough to pro¬ 
vide _a robust exclus ion of the parameters inferred by 


Bulbul et al. (2014a), a Micro-X observation of the GC 
or a region near the GC could provide a significantly 
stronger constraint because of its better energy resolu¬ 
tion and larger exposure, and because of the larger sig¬ 
nal strength in the GC region. These limits obviously 
depend on the structure of the dark matter halo, with 
cored profiles producing stronger constraints than NFW- 
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Figure 12. Constraints on decaying sterile neutrino dark matter, assuming that sterile neutrinos comprise all of the DM in the NFW 
profile of|Nesti & Salucci| (20! d|). Limits include the XQC observation analyzed in this work (black)-, Micro-X median expectation from an 
observation of the L.L (1, doWed blue), an observation below the plane of the galaxy in the direction of l = 0 °,b = —32° (2, sho rt dashed 
blue), and an observation within the XQC field in the direction of l = 162°, b = 7° (3, long dashed blue h con straints from M31 (Horiuchij 
let al.l20 14| (sha ded ora nge)-, and constraints from the previous analysis of XQC data by|Boyarsky et al.|l|2007|l (gray). The putative signal 
offfjiilbiil ot al.j (2014a;) is also shown (red point). 


like profiles. 

These wide-FOV rocket observations are complemen¬ 
tary to narrow-FOV observations of dwarfs, galaxies, and 
clusters because they directly address whether an uniden¬ 
tified line is present as an all-sky signal in the MW. 
This confirmation would be crucial for distinguishing an 
atomic interpretation from an exotic DM one, and for 
establishing the signal scaling as a function of the inte¬ 
grated DM density. 

6. CONCLUSION 

Microcalorimeters onboard sounding rockets have the 
ability to place competitive bounds on keV sterile neu¬ 
trinos or other dark matter models whose flux scales 
linearly with dark matter density. We have analyzed 
a subset of the data acquired during the 5th flight of 
the XQC payload corresponding to an effective exposure 
of 106 s on 24 pixels and placed a upper limit on keV 
sterile neutrinos between 4-10 keV which demonstrates 
the prospect for future observations with this type of in¬ 
strument. A study of future observations in and around 
the Milky Way galactic center with the Micro-X pay- 
load shows that it will have sensitivity to new parameter 
space in the (m s , sin 2 20) sterile neutrino plane. Op¬ 
timizations of the pointing direction, increased field of 
view and energy resolution, and repeated observations 
will all increase the sensitivity of this technique in the 
future. 
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APPENDIX 


A. STATISTICAL MODEL FOR XQC 

We use the method of unbinned extended maximum 
likelihood to fit the XQC data. The unbinned method 
produces equivalent results to binned y 2 fits in the limit 
of high statistics, but produces more reliable fits in low 
statistics set tings where many bins would have few or 
zero events (James 20061. The likelihood function for 
our background model has the form 


£(S -,{Ei}) = 


e~^n N 

N\ 


N 

n 


i= 1 lk=l 


y;—^(s fc ;^) 


, (Al) 


where the product is taken over the N total events in the 
observation, and the sum is taken over each of 7 compo¬ 
nents of the background model. The values fik are the 
estimated number of events in each of the background 
components and fi = The probability density 

functions (PDF) for each component are the Pk( S*,; Ei), 
which are functions of energy and depend on the vec¬ 
tor of parameters S*. (with S = 1% SC. The functional 
forms for each PDF are listed in Table [Al] and the corre¬ 
sponding parameters are in Table fA2| Background PDFs 
which have a fixed template shape, such as the cosmic 
rays, do not have any parameters, so S /. is an empty 
set. The Poisson term (first) is the extended likelihood 
term, which constrains the total expected event rate by 
the number of observed counts. 

The compon ents of the background PDF are summa¬ 
rized in Table |A1[ an d the key model parameters are 
contained in Table I A2l The calibration lines from inter¬ 
actions in the HgTe absorber are modeled using a gaus- 
sian kernel density estimate (KDE) based on calibration 
data taken in a lab after launch. An alternate model 
using Voigt profiles for each of the calibration lines also 
produces a reasonable fit, but the KDE-based model was 
chosen because it agrees better in the low-energy tails 
of the calibration lines. A similar KDE does not accu¬ 
rately describe the corresponding lines in the substrate, 
around 2.80 keV and 3.0 keV. These are more reliably 
modeled by gaussian PDFs with fitted means and a com¬ 
mon fitted width. Events that interact in the substrate 
have suppressed energy because of charge trapping ef¬ 
fects that depend on the neutralization state of the Si, 
so differences between flight and calibration data are not 
surprising. Lastly, the photoelectron produced by an X- 
ray interaction in the HgTe absorber escapes the active 
detector volume in about ~ 5% of events. We model 
this by a flat distribution extending from zero energy to 
the Kq line, whose normalization is fixed to 5% of the 
number of events in the Ka and K/3 absorber lines. 

We model the X-ray continuum with two power laws. 
One describ es the flux fr om the Crab, using canonical 
parameters (Mori et al. 2004). The other describes the 
cosmic X-ray background, using parameters f rom the 
Chandra deep field measurement of |Hickox fc Markevitch] 
(2006). Before fitting, the power laws must be weighted 
by the efficiency of both the HgTe absorber and the Si 
substrate. Since the events in the Si substrate appear 
below their true energy, this component must be shifted 
to lower energies by a similar fractional energy loss as the 



i i I i i i I i i i I i i i I i i i l , 1 1 ‘cr-rMVr-l 

0 2 4 6 8 10 12 14 16 18 20 


Energy [keV] 


Figure Al. Probability density function of simulated energy 
spectrum in the X QC (blue) and Micro-X (dashed green) X-ray 
absorbers. Geant4 (jAgostinelli et al.||2003( is used t o simulate a 
powe r law distribution ot primary cosmic ray protons (|Papini et al.| 
|1996|) impinging isotropically on the two absorbers. Both distribu¬ 
tions are typical for minimum-ionizing particles. The mean energy 
deposited in XQC is larger than in Micro-X because of the addi¬ 
tional 15 /am Si substrate of the HgTe absorber present in XQC 
but not in Micro-X. 


calibration lines. The resulting PDF for the continuum 
is given by 


Ppower(E) = e HgTe (E)E 


~ a + e subst ( E / k ) ( ^ 


(A2) 

where k = E^ st - / 3.31 keV is an estimate of the frac¬ 
tional energy loss of the K a calibration line. Note that 
this parameterization implicitly assumes that the charge 
trapping process in the substrate is energy-independent. 
Since the Crab lies 19.5° off the observation axis, the ge¬ 
ometrical acceptance of the detector to X-rays from the 
Crab is 94.2% of the effective area for on-axis events. 
Because of a similar geometrical effect, the acceptance 
of the diffuse X-rays is 92.7% of the effective area for 
on-axis events. 

The background component due to cosmic rays is only 
about 5 events in the 2-5 keV window for the full XQC 
exposure of 2551 pixel s. We obtain the spectral shape of 
cosmic rays by simulating proton s with a typical power 
law spectrum (Papini et al.| 1996) impinging on the XQC 
absorber and substrate. Since protons are minimum ion¬ 
izing particles, this is approximately a Landau dist ribu - 
tion with a peak around 7 keV, as shown in Figure [Al] 
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k 

PDF component 

Functional form 

1 

potassium Kq, K/3 cal. events in HgTe 

KDE based on pre-flight calibration 

2 

potassium Ka cal. events in Si substrate 

gaussian, with fitted mean and width 

3 

potassium K/3 cal. events in Si substrate 

gaussian, with fitted mean and width 

4 

photoelectrons from cal. source escaping from absorber 

uniform from 0 keV to potassium Ka energy 

5 

Crab nebula 

power law 

6 

cosmic X-ray background (CXB) 

power law 

7 

cosmic rays 

spectrum derived from Geant4 simulation 
of protons with power law above 1 GeV (cc = 2.7) 


Table Al 

Components of XQC background PDF in equation ED- 


k 

Parameter 

Value 

i 

number of Ka, K/3 cal. events in HgTe 

1281±36 

2 

number of Ka cal. events in Si substrate 

722 ± 29 

2 

mean measured energy of Ka cal. line in Si substrate 

2.784 ±0.006 keV 

2 & 3 

width of measured Ka, K/3 cal. lines in Si substrate 

0.061 ± 0.002 keV 

3 

number of K/3 cal. events in Si substrate 

221 ± 18 

3 

mean measured energy of K/3 cal. events in Si substrate 

3.020 ± 0.006 keV 

4 

number of cal. source events with escaping electrons in fit range 

25 (fixed 5% escape fraction) 

5 

Crab spectral index (Mori et al.|2004 i 

Crab flux at 1 keV (Mori et al. 2004) 
number of Crab events 

2.1 (fixed) 

9.7 ± 0.5 photons cm -2 s -1 keV -1 (fixed) 

155 (fixed) 

6 

CXB spectral index (Hickox Sz Markevitch)2006) 
constraint on CXB flux at 1 keV ([Hickox Sz Markevitch|2006 i 
number of CXB events 

1.4 (fixed) 

10.9 ± 1.3 photons cm -2 s _1 keV -1 sr -1 (fixed) 
369 (fixed) 

7 

number of cosmic rays in fit range 

5.2 (fixed) in 2-5 keV range 

— 

exposure time 

2551 pixel s (106 s on 24 pixels) 


Table A2 

Parameters of the XQC background PDF in equation |aTJ. 
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